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Abstract. The probability density functions (PDFs) for energy dissipation rates, 
created from time-series data of grid turbulence in a wind tunnel, are analyzed in a high 
precision by the theoretical formulae for PDFs within multifractal PDF theory which 
is constructed under the assumption that there are two main elements constituting 
fully developed turbulence, i.e., coherent and incoherent elements. The tail part of 
PDF, representing intermittent coherent motion, is determined by Tsallis-type PDF 
for singularity exponents essentially with one parameter with the help of new scaling 
relation whose validity is checked for the case of the grid turbulence. For the central 
part PDF representing both contributions from the coherent motion and the fluctuating 
incoherent motion surrounding the former, we introduced a trial function specified 
by three adjustable parameters which amazingly represent scaling behaviors in much 
wider area not restricted to the inertial range. From the investigation of the difference 
between two difference formulae approximating velocity time-derivative, it is revealed 
that the connection point between the central and tail parts of PDF extracted by 
theoretical analyses of PDFs is actually the boundary of the two kinds of instabilities 
associated respectively with coherent and incoherent elements. 
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1. Introduction 

There are several keystone works (Mandelbrot 1974, Parisi and Frisch 1985, Benzi et al 
1984, Halsey et al 1986, Meneveau and Sreenivasan 1987, Nelkin 1990, Hosokawa 1991, 
Benzi et al 1991, She and Leveque 1994, Dubrulle 1994, She Z-S and Waymire 1995, 
Arimitsu T and N 2000a, 2000b, 2001, 2002, Arimitsu N and T 2002, Biferale et al 2004, 
Chevillard et al 2006) providing the multifractal aspects for fully developed turbulence. 
Only a few works (Benzi et al 1991, Arimitsu T and N 2001, 2002, Arimitsu N and T 
2002, Biferale et al 2004, Chevillard et al 2006) analyze the probability density functions 
(PDFs) for physical quantities representing intermittent character. The other works deal 
with only the scaling property of the system, e.g., comparison of the scaling exponents 
of velocity structure function. Among the researches analyzing PDFs, multifractal 
probability density function theory (MPDFT) (Arimitsu T and N 2001, 2002, 2011, 
Arimitsu N and T 2002, 2011) provides the most precise analysis of the fat-tail PDFs. 
MPDFT is a statistical mechanical ensemble theory constructed by the authors (T.A. 
and N.A.) in order to analyze intermittent phenomena providing fat-tail PDFs. 

To extract the intermittent character of the fully developed turbulence, it is 
necessary to have information of self-similar hierarchical structure of the system. This is 
realized by producing a series of PDFs for responsible singular quantities with different 
lengths 

C = £o5-", 5>1 (n = 0,l,2,---) (1) 

that characterize the regions in which the physical quantities are coarse-grained. The 
value for 5 is chosen freely by observers. Let us assume that the self-similar structure of 
fully developed turbulence is such that the choice of 8 should not affect the theoretical 
estimation of the values for the fundamental quantities characterizing the turbulent 
system under consideration. A&A model within the framework of MPDFT itself tells 
us that this requirement is satisfied if the scaling relation has the form (Arimitsu T and 
N 2011, Arimitsu N and T 2011) 

In 2/(1 -q)\n5= l/a_ - l/«+. (2) 

Here, q is the index associated with the Renyi entropy (Renyi 1961) or with the Havrda- 
Charvat and Tsallis (HCT) entropy (Havrda and Charvat 1967, Tsallis 1988); a± are 
zeros of the multifractal spectrum f{a) (see below in section [2]). The multifractal 
spectrum is uniquely related to the PDF for a (see (jl]) below). The PDF of a is related 
to the tail part of PDFs within MPDFT for those quantities revealing intermittent 
behavior whose singularity exponents can have values a < 1, e.g., the energy dissipation 
rates, through the variable transformation between a and the physical quantities (see ([3]) 
below for the case of the energy dissipation rates e„). With the new scaling relation (j2]), 
observables have come to depend on the parameter 5 only in the combination (1 — q') In 5. 
The difference in 5 is absorbed in the entropy index q. |§| 

§ Since almost all the PDFs that had been provided previously were for the case where (5 = 2, it 
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In the preceding papers, we analyzed PDFs for energy transfer rates (Arimitsu T 
and N 2011) and PDFs for energy dissipation rates (Arimitsu N and T 2011), which 
are given in figure 11 of Aoyama et al (2005), with the help of the new scaling relation, 
and checked the independence of the PDFs from 5. It was found that the adjustable 
parameters for the central part PDF provide us with 5-independent scaling behaviors 
as functions of r/r/, and that the scaling properties are satisfied in much wider region 
not restricted to inside of the inertial range. However, the number of data points used 
in drawing the PDFs in figure 11 of Aoyama et al (2005) is not enough, especially, for 
the precise analyses of the central part of the PDFs performed in Arimitsu T and N 
(2011) and Arimitsu N and T (2011). Therefore, we will perform, as one of the aims 
of the present paper, the same analyses, which were done for DNS, with the help of 
PDFs created from wind tunnel turbulence with a higher enough resolution in order 
to make sure if the characteristics discovered previously with rather poor resolution at 
the central part are correct or not. Since we have the raw time-series data taken from 
wind tunnel turbulence, we can create PDFs for energy dissipation rates with enough 
resolutions fit to our needs. 

In this paper, we analyze the PDFs for energy dissipation rates extracted out 
from the time series of the velocity field of a fully developed turbulence which were 
observed by one of the authors (H.M.) in his experiment conducted in a wind tunnel 
(Mouri et al 2008). In section O we present the formulae of theoretical PDFs within 
A&A model which are necessary in the following sections for the analyses of PDFs 
obtained from the experimental turbulence. In section [3l we analyze the observed PDFs 
for energy dissipation rates in a high precision with the theoretical PDF within A&A 
model of MPDFT, and verify the proposed assumption related to the magnification 5. 
In section |U in order to see what information we can extract out from the time-series 
data, we compare two different PDFs for energy dissipation rates created from the time 
series data with different approximation for temporal derivative. We may learn from this 
how to treat the central part of PDFs to derive the information of incoherent fluctuating 
motion around the coherent turbulent motion. Summary and discussion are provided 
in section [51 

2. Singularity exponent and PDFs for energy dissipation rates 

MPDFT is constructed under the assumption, following Parisi and Frisch (1985), that 
for high Reynolds number the singularities distribute themselves in a multifractal way 
in real physical space. The singularities stem from the invariance of the Navier-Stokes 
(N-S) equation for an incompressible fluid under the scale transformation x — )■ x' = Ax, 

has been possible to analyze PDFs (Arimitsu T and N 2001, 2002, Arimitsu N and T 2002) with the 
scaling relation 1/(1 — (?) = l/a_ — l/a+ proposed by Costa et al (1998) and Lyra and Tsallis (1998) in 
connection with the 2°° periodic orbit. The orbit having the marginal instability of zero Liapunov 
exponent appears at the threshold to chaos via a period-doubling bifurcation in one-dimensional 
dissipative maps. 
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accompanied by the scale changes u ^ u' = A"/^m in velocity field, t t' = A^~"/^t in 
time and p p' = X^^^^p in pressure with an arbitrary real number a, in the limit of 
large Reynolds number, i.e., the contribution from the dissipation term in N-S equation, 
which is proportional to the kinematic viscosity z/, is negligibly small compared with the 
convection term. In treating an actual turbulent system, however, the value u is fixed 
to a finite value unique to the material of fiuid prepared for an experiment. We should 
keep in mind that the dissipation term can become effective depending on the region 
under consideration, since the term breaking the invariance does exist, i.e., non-zero 
(see the discussion in the following). 

The invariance under the scale transformation leads to the scaling property 

Sn/e = {ejior-' (3) 

for the energy dissipation rate averaged in the regions with diameter in- Here, we 
put Eq = e whose value is assumed to be constant. The energy dissipation rate becomes 
singular for a < 1, i.e., hm„_j.oo£^n = lim„_s>ooC~^ ~^ The degree of singularity is 
specified by the singularity exponent a (Parisi and Frisch 1985). 

Let us consider a to be a stochastic variable whose PDF P^"'\a) is given by the 
Renyi or HCT type function (Arimitsu T and N 2000a, 2000b, 2001, 2002, 2011, Arimitsu 
N and T 2002, 2011): 

P(")(a) oc [1 - (« - ao)V(Zi«)T^^'"'^ (4) 

1 /2 

with Aa = [2X/(1 — q)\n6] ' . The domain of a is < a < Omax with amin and 
Omax being given by 

Q^min/max = ao =F ^a- Q is the entropy indexIJi] From (j4]), we have 
for n ^ 1 the expression of the multifractal spectrum 

/(a) = 1 + In [1 - (a - aof/{Aaf] /{l-q) In 6. (5) 

The independence of /(«) from n is interpreted as a manifestation of the existence of self- 
similar hierarchical structure responsible for the intermittent fiuid motion of turbulence. 

The three parameters Uq, X and q appeared in P^^\a) are determined as the 
functions of the intermittency exponent /i with the help of the three conditions. One 
is the energy conservation law = e. Another is the definition of the intermittency 
exponent /x, i.e., ((£:„/e)^) = (^n/^o)"^- The last condition is the scaling relation ([2]) 
with aj- being the solution of f{a±) = 0, which is a generalization of the one introduced 
by Tsallis and his coworkers (Costa et al 1998, Lyra and Tsallis 1998) to which ([2]) 
reduces when 6 = 2. Here, the average (■ ■ ■) is taken with P^'^\a). The parameter q 
is determined, altogether with and X, as a function of /i only in the combination 

II The function ^ is the MaxEnt PDF derived from the Renyi entropy or from the HCT entropy 
with two constraints, one is the normahzation condition and the other is a fixed q-variance (TsalUs 
1988). This choice of PDF is also quite natural since the Renyi entropy and the HCT entropy are 
directly related to the generalized dimension (Hentschel and Procaccia 1983) describing those systems 
containing multifractal structures (Grassberger 1983). Note that for the HCT entropy the relation is 
given with the help of the g-exponential (Tsallis 2001) which is a function satisfying a scaling invariance 
(Suyari and Wada 2006) and reduces to the ordinary exponential for q ^ 1. 
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(1 — g) ln(5. The difference in S is absorbed into the entropy index q, therefore changing 
the zooming rate 6 may resuh in picking up the different hierarchy, containing the 
entropy specified by the index q, out of self-similar structure of turbulence. As the 
parameters are dependent on q only in the combination (1 — q)\nS, we are naturally 
led to the replacement of n in the expression of P^^^{a) in (jlj) with n = n/lia.5. If h 
does not depend on 6, P^''\a) becomes also independent of Note that, with the 
new number n, introduced in ([1]) reduces to 

L = ^oe-^. (6) 

MPDFT provides us with a systematic framework to make a connection between the 
PDF P*^"-'(a) of the singularity exponent a and the PDF of the observed quantity such 
as the energy dissipation rate representing intermittent singular behavior in its time- 
evolution. The element of fluid motion specified by the singularity exponent satisfying 
a < 1 takes care of the intermittent large (singular) spikes observed in the time-evolution 
of energy dissipation rate, and contributes to the tail part of PDF for energy dissipation 
rates (see figured] (a) and (b)). This element is directly related to a coherent hierarchical 
structure such as the multi-scale Cantor set characterized by the multifractal spectrum 
/(a). Therefore, the fluid motion controlled by this element is referred to as a coherent 
motion. There is another element of fluid motion due to the symmetry breaking term, 
i.e., the dissipation term in N-S equation, which produces fluctuation of fluid surrounding 
the coherent turbulent motion. This element contributes mainly to the central part of 
PDF (see figured] (a) and (b)). The fluid motion provided by this element is referred to 
as an in-coherent motion. Note that the central part of the PDF is constituted of two 
elements, i.e., the in-coherent and coherent motions. 

Based on the above consideration, we assume that the probability I\.^^\en)den can 
be, generally, divided into two parts as 

n^^\en)den = ng (£„)&„ + Ant\en)den (7) 

(see figure 1 (a) and (b)). The first term describes the coherent motion, i.e., the 
contribution from the abnormal part of the physical quantity e„ due to the fact that 
its singularities distribute themselves multifractal way in real space. This is the part 
representing a coherent turbulent motion given in the limit — )■ but is not equal 
to zero (z/ ^ 0). The second term represents the contribution from the incoherent 
fluctuating motion. The normalization of PDF is specified by derJl^^\sn) = 1- 
We assume that the coherent contribution is given by (Arimitsu T and N 2001) 
^3^s(^n)den = T[^^^P'^^\a)da with the variable transformation (|3]). For the expression 
of n^"c], see Arimitsu N and T (2011). 

% The introduction of h is intimately related to the infinitely divisible process (Dubrulle 1994, She 
and Waymire 1995). It is confirmed by the observation in the preset paper that n is independent of 5 
and has values of 0(1) (see table [5]). Then, taking the limit 5-^1+ with a fixed value of n, one has 
an infinitely divisible distribution. A detailed investigation of A&A model from this view point will be 
given elsewhere in the near future. 
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Figure 1. Two kinds of divisions of PDF Ilg ^(e„). One into Ilg g(e„) and Z\n3 '(£„) 
are given on (a) linear and (b) log scale in the vertical axes. The other into fl'i^li^n) 
and n3"''i(Cn) is given in (c) on log scale. The open circles represent an experimental 
PDF for energy dissipation rates. The contribution of Z\n3"''(e„) to the tail part 
ns'til^i) is negligibly small. 



Let us introduce another division of the PDF (see figure 1 (c)), i.e., 

ni"Hen) = n£(en) + n£l(a), (8) 

where n3"''(^„) is introduced by the relation fl^^\c,n)dC,n = 1^3"^ with the variable 

transformation ^„ = e„/ ((e^))^^ where the cumulant average ((■ ■ ■))c is taken with the 
PDF 113"^ The two parts of the PDF, the tail part Il^^^i{^n) and the central part 
^t^cii^n), are connected at = C under the conditions that they have the common 
value and the common log-slope there. Note that ^* is related to e* by ^* = e* / {{£n))c^'^ 
and to a* by (IH]). The value of a* is determined for each PDF as an adjusting parameter 
in the analysis of PDFs obtained by ordinary or numerical experiments. 

When one creates a PDF from the time-evolution data for microscopic energy 
dissipation rate, he puts each realization into an appropriate bin according to the value 
En which is obtained by averaging the microscopic energy dissipation rates in each time 
interval corresponding to the length in- For larger En values belonging to the tail part 
domain of the PDF, most of the realizations in a bin at the interval ~ + come 
from the time interval containing at least one intermittently large spike (singular spike) 
of microscopic energy dissipation rates. The bin may have negligibly small proportion of 
the number of realizations coming from those intervals with only fluctuations compared 
to the number of realizations with at least one singular spike. On the other hand, for 
smaller En values belonging to the central part PDF domain, the number of realizations 
coming from the time intervals containing singular spikes with smaller heights is about 
the same order as the number of realizations from the time intervals containing only 
fluctuations, since the height of singular spikes contributing to this bin must be about 
the same height as fluctuations. 

Under the above interpretation, it may be reasonable to assume that, for the 
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tail part of PDF Il^^■^{^n), the contribution from the first term Il.^a^{en) in ([7]) to the 
intermittent rare events dominates, and the contribution from the second term Zillg"'' (5„) 
to the events is neghgible, i.e., 

for ^* < C,n- For < ^„ < as there is no theory for the central part of PDF n3"'*j,(^„) 
at present, we put 

n£(en)rfen = n(")e-[^3(c.)-.3(«rji (e„/c) d^n (10) 

with 113"^ = fl^^l^y\f"{ao)\/2^T\ ln(£„/£o)|/^n and a trial function of the Tsallis-type 

X {1 - (1 - g') [9 + /'(a*)] [(a/C)'"^ - 1] M'^^'-''^ (11) 

containing minimal number of adjustable parameters, i.e., q', 9 and w^. The parameter 
is adjusted by the property of the experimental PDFs around the connection point; 
q' is the entropy index different from g in (jl]); 6' is determined by the property of PDF 
near ^„ = 0. For the expression of see Arimitsu N and T (2011). The contribution 
to n3"jj,(^„) comes both from coherent and incoherent motions (see figure [1]). 

The reason why we chose the trial function (fTT!) for the central part PDF is because 
it is a natural generalization of the x-square distribution function for the variable 
Un = (Cn/C)'"^- The observed value of q' is in the range 1.03 < q' < 1.09 (see table [2]). 
Note that in the limit g' — ?• 1 the trial function reduces to the x-square distribution 
function for The quantity {9 + ws — l)/ws provides us with an estimate for the 
number of independent degrees of freedom for the dynamics contributing to the central 
part of PDF. 

3. Verification of assumptions through the analyses of experimental PDFs 
3.1. Experimental setup and extraction of PDFs 

By means of the theoretical formula within MPDFT summarized in the last section, we 
analyze PDFs of energy dissipation rates created from the time series data (Mouri 
et al 2008) for the turbulence produced by a grid in a wind tunnel (see table [1]). 
Measurements are performed by a hot-wire anemometer with a crossed-wire probe 
placed on the centerline of the tunnel at 4 m downstream from the grid. It is expected 
that turbulence around the probe is homogeneous in both stream-wise and span-wise 
directions, as the cross-section 16 cm x 16 cm of each open square surrounded by the 
rods constituting the grid is small enough compared with the cross-section 3 m x 2 
m of the wind tunnel. The cross section of a rod is 4 cm x 4 cm. We also expect 
that the turbulence around the probe is isotropic even for larger scales since the values 
of RMS one-point velocity fluctuations for span-wise and stream-wise components are 
almost equal (see table 1). There are still possible pitfalls about the assumption of 
isotropy (Biferale and Procaccia 2005) because of the difference of values between the 
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Table 1. Parameters of the grid turbulence in a wind tunnel (Mouri et al 2008). The 
inertial range is determined as the region where the second moment of the velocity 
differences for longitudinal component scales with the exponent 2/3 with respect to the 
distance between the positions of two velocities used to derive the velocity difference. 



Quantity 



Value 



Microscale Reynolds number Rca 
Kolmogorov length 77 
Kinematic viscosity u 
Mean velocity of downstream wind U 
Mean energy dissipation rate ((e)) = 15u{{{dv/dx)'^)) /2 
Correlation length of longitudinal velocity 
Inertial range 
RMS of span- wise velocity fluctuations ((w^))^^^ 
RMS of stream- wise velocity fluctuations ((w^))^^^ 
Sampling interval At 
Number of data points 



1.42 X 10"^ m^ sec 



50 < r/T] < 150 
1.06 m sec~^ 



1.10 m sec 
1.43 X 10~5 
4 X 10* 



21.16 m sec ^ 
7.98 m^ sec-3 



0.138 mm 



17.9 cm 



409 



-1 



sec 



-1 



averaged energy dissipation rates estimated with the span-wise velocity component v, 
i.e., 15z/((((9f /9x)^))/2 = ((e)) = 7.98 m^ sec^^, and the one estimated with the stream- 
wise velocity fluctuation u, i.e., 15z/(((9m/9x)^)) = 8.58 m^ sec~^ (Mouri et al 2008). 
However, as the difference is less than 10 %, we expect that anisotropy, even if it exists, 
may not affect the following analyses seriously. 

Assuming isotropy of the grid turbulence, we adopted the surrogate 
lhv{dv / dxY /2 = 15u(dv/dty /2U'^ for the energy dissipation rate (Cleve et al 2003, 
Mouri et al 2008) with the mean velocity U of downstream wind (see table [1]) where 
X-axis is chosen to the direction of the mean flow in a wind tunnel and v is the span-wise 
velocity component. Here, we used Taylor's frozen hypothesis in replacing the variable 
from time t to space x (see Mouri et al (2008) for detail). For the estimation of dv/dt, 
we use here the difference formula 



where At is the sampling interval observing velocity (see tabled]). With this formula, we 
can have a better estimate of the velocity time derivative by means of A^^h/At without 
contamination up to the term of 0{AtY. We represent the /oca/ energy dissipation rates 
derived from ( IT2l) by the symbol e, i.e., e = 15z/(Z\*^^^f /Z\t)^/2f/^. 

In creating the experimental PDFs for energy dissipation rates, 4 x 10^ data points 
are put into 2 x lO'^ bins along the ^„ axis. We discarded those bins containing the 
number of data points less than 25. Note that the average number of data points per 
bin is 2 X 10^. In drawing the created PDFs for energy dissipation rates, not all the 
bins but every 10^ bins are plotted for better visibility. The experimental PDF in the 
region near the right-most end points are scattered because of the lower statistics due to 
smaller number of data points in the bins located there (see figure [2] (a) and figure H] (a)). 



dv/dt ~ A^^^v/At = {8 [v{t + At) - vit - At)\ 

-\v{t + 2At) - v{t - 2At)]} /UAt 



(12) 
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Figure 2. PDFs of energy dissipation rates for (5 = 3 on (a) log and (b) linear scale 
in the vertical axes. For better visibility, each PDF is shifted by —2 unit along the 
vertical axis in (a) and by —0.4 unit along the vertical axis in (b). Closed circles are the 
experimental PDFs for r/rj — 21.9, 65.7, 197 and 591 from the smallest value (top) to 
the largest value (bottom) where r corresponds to £„. Solid lines represent the curves 
given by the present theory with parameters listed in table [2] (a) . 



3.2. Analyses of experimental PDFs 

The experimental PDF is analyzed with the help of the theoretical formula for PDF 
by the following procedure: (i) Pick up three experimental PDFs with consecutive r 
values, say, ri, r2 = ri5 and rs = ri5^. (ii) With a trial /x value, analyze each of the 
three experimental PDFs to find out tentative but the best values g', W3, 6*, a* and 
Ui = ln(rj/£o)/ln5 [i = 1, 2, 3) for the theoretical PDF. (iii) Check if the differences 
77-3 — n2 and n2 — ni are close to 1 or not. (iv) If not, change /i value, and repeat 
the processes (ii) and (iii) until one arrives at the set of best fit parameters under the 
condition ns — n2 = n2 — ni ~ 1 within a settled accuracy, (v) With thus determined 
common /i value, determine the best fit values q', w^, 9 and a* for each of other PDFs 
which are not picked out for the above processes (i) to (iv). One notices that rij— nj_i ~ 1 
are satisfied automatically for every PDFs created from the experiment. 

The PDFs of energy dissipation rates are analyzed in figure |2] for the magnification 
(5 = 3 on (a) log and (b) linear scale in the vertical axes. For better visibility, each PDF 
is shifted by appropriate unit along the vertical axis. Closed circles are the experimental 
data points for PDFs for the cases r /t] = 21.9, 65.7, 197 and 591 from the smallest value 
(top) to the largest value (bottom) where r corresponds to Solid lines represent the 
theoretical PDFs given by ([H]) with (jS]) and ( fTOj) . The parameters necessary for the 
theoretical PDF of A&A model are determined as {1 — q)\n5 = 0.393, = 1.15 and 
X = 0.310, which turn out to be independent of S. Other parameters are listed in 
table E] (a) and table [3] (a). We performed the same analyses for other magnifications, 
(5 = 2 and 5, and found that the extracted value fi = 0.260 for the intermittency 
exponent is common to three cases in which PDFs are created with the different values 
of magnification, i.e., 5 = 2, 3, 5. It means that, within the analysis of the energy 
dissipation rates, the turbulent system under consideration is characterized by a unique 
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Table 2. Parameters of PDFs created by (a) the formula (|12p and (b) the formula 
([T3| . For both cases, ^ = 0.260 ((1 - q)\nS = 0.393, = 1.15, X = 0.310) giving 
q = 0.642. 
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2.20 
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0.750 
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0.60 


0.416 


1.09 


1.25 


5.80 


0.580 


0.637 


1.09 


1.24 


6.20 



Table 3. Connection points between the central and the tail part PDFs created by 
(a) the formula ^ and (b) the formula ((e)) = 7.98 m^ sec-^. 
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5.24 
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value as it should be. 




0.5 







(b) 


- 








10- 










\\* ■ 






























5- 




®^ 


A 










0- 



10 



10' 



10- 




Figure 3. r/rj (= inM) dependence of (a) (b) a* and (c) 6. In (a), the data points 
extracted by the present analysis via are plotted by closed circles for ^ = 2, by 
closed squares for 5 = 3, by closed triangles for 5 = 5, whereas those extracted by the 
DNS analysis (Arimitsu N and T 2011) are plotted by symbols nabla for 5 — 2^/^, 
by times for 5 = 2^/^, by pluses for 5 = 2. The empirical formula for the present 
(DNS) analysis is given by n = -2.39 logio(r-/r7) + 7.31 (n = -2.33 logio(?'/?7) + 8.74). 
In (b) and (c), the data points extracted by the present analyses via (fT2|) (via 
([T3|) ) are plotted by closed (open) circles for (5 = 2, by closed (open) squares for 
5 = 3, by closed (open) triangles for (5 = 5. Solid (dashed) lines are the empirical 
formulae (b) a* = -0.326 logio(r/?7) + 1.05 (a* = -0.285 logio(7-/?7) + 0.966) and (c) 
e = 1.831ogio(?'/?7) + 0.460 {9 = 1.32 logio(r/ry) + 2.70). The empirical formulae are 
obtained by using all the data points for different values of 5. The inertial range for the 
present (DNS) analysis is the region between the vertical dash-dotted (dotted) lines. 



The dependence of n on r/r] (= £„/ rj) for the present analysis with the series of PDFs 
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derived through (fT^ is given in figure [3] (a) by closed circles for 5 = 2, by closed squares 
for S = 3 and by closed triangles for 6 = 5. The empirical formula for n obtained by 
making use of all the data points for 5 = 2,3 and 5 with the method of least squares has 
the expression n = — 1.031n(r/?7) + 7.31 which is drawn by a solid line (lower line in the 
figure). This proves the correctness of the assumption that the fundamental quantities 
of turbulence are independent of 6. We also include in the figure, for comparison, the 
data points of n for 4096^ DNS taken from figure 4 in Arimitsu N and T (2011) and the 
empirical formula n = —1.01 ln(r/?7)+8.74 (upper solid line) derived with the data points 
by the method of least squares. For the DNS, /i = 0.345 (Arimitsu N and T 2011). How 
much n data points are scattered from the empirical formula (see figure |3] (a)) and also 
from the theoretical formula ([6]) with = r provides us with a measure how much we 
perform appropriate extraction of parameters. The data points of n for the turbulence 
in the wind tunnel are scattered more compared with those for the turbulence in 4096^ 
DNS, as the time-series raw data from wind tunnel include indispensable measurement 
errors associated with readout processes, e.g., mainly, electrical noises. 

The r/r] (= in/v) dependences of a* and 9 are given, respectively, in figure [3] 
(b) and (c) by closed circles for 5 = 2, by closed squares for 5 = 3 and by closed 
triangles for 6 = 5, which are extracted from the series of PDFs derived through (fT2|) . 
The solid line in each figure, (b) and (c), represents an empirical formula obtained 
from all the data points for 6 = 2, 3 and 5 by the method of least squares. These 
figures prove again the correctness of the assumption that the fundamental quantities 
of turbulence are independent of 6. The value of q' is found to be about q' = 1.05 
(see table [2] (a)). We found that is also independent of 6 and has a common line 
logio ^3 = 0.372 logiQ^r /rj) + log^g 0.112. Note that the empirical formulae for n, a* and 
6 is effective only for the region r/rj >2 since 6 should satisfy ^ > 1 (see figure [3] (c)). 
We observe that the parameters q', 6, for the central part PDF and the connection 
point a* have scaling behaviors in much wider region not restricted to inside of the 
inertial range. 

4. Comparison of PDFs produced with full and less contaminations 

In this section, we analyze the PDFs for the energy dissipation rates derived from the 
time-series data with the difference formula 

dv/dt ~ Z\(°)i;/Z\t = [v{t + At) - v{t)] /At (13) 

in order to study what difference comes out compared with the PDFs analyzed in 
section [3] which is derived by means of the difference formula (JT2l) . Note that the 
formula f|T3l) estimates the values of velocity time derivative with A^^^v/At which 
may contain full contamination, i.e., from the 1st order term with respect to At. 
We introduce the symbol for the local energy dissipation rates derived from (fT3|) . 
i.e., e^^^ = 15h'{A^^^v / At)'^ /2U'^^ In creating the experimental PDFs for the energy 

+ We observe that ((e'-'^-')) = 8.08 m'^ sec^"'^ which is larger than ((e)). 
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Figure 4. Comparison of PDFs for energy dissipation rates Ilg (e/((e))) and 
n3"^(e'°V((^))) created, respectively, with the formulae (IT^ and ([T^ . In (a) and 
(d)-(f), closed (open) circles and full (dashed) lines are, respectively, the experimental 
and theoretical PDFs for n^"^(e/((£))) {n'^\e^°y {{e)))) with ^ = 0.260. PDFs in (a) 
represent for the cases r/rj = 6.57 (top), 21.9 (middle) and 65.7 (bottom), shifted 
by —2 unit along the vertical axis for better visibility. The magnification of the 
central part PDFs for each r (= in) is given in (d)~(f). The relative difference 
A„ = [n(") (£(")/((£))) - (£/ ((e)) )]/4") (£/((£))) is given for (b) r/ry = 6.57 and (c) 
65.7, in which closed circles (full lines) are experimental (theoretical) A„. The error 
bar is the standard deviation of 100 hidden (not appeared in the figures) data points 
for A„ which locate between the adjacent data points for A„ appeared in the figures. 
The shown error bars are thinned out. Note that (a) and (b)-(f) are, respectively, 
drawn on log and linear scale in the vertical axes. 

dissipation rates, we took the same procedure as used in section |3l 

We compare, in figures H] (a) and (d)-(f), the PDFs of energy dissipation rates 
Il^^\e / {{e))) and Il^^\e^^'^ / {{e))) which are created, respectively, with the help of formulae 
( IT2I) and ( IT3I) . Note that the arguments for every PDFs are scaled by ((e)) which 
does not depend on r (= in)- In figure H] (a) each PDF is displayed on log scale in 
vertical axis for the cases r/r] = 6.57 (top), 21.9 (middle) and 65.7 (bottom), which 
are shifted by —2 unit along the vertical axis for better visibility. The magnification of 
their central part PDFs are displayed in figures ll](d) r/rj = 6.57, (e) 21.9 and (f) 65.7 
on linear scale in vertical axis. The closed (open) circles and the full (dashed) lines are, 
respectively, the experimental and theoretical PDFs for Il^^\e / {{e))) {Il^^\e''^'' / {{e)))) 
with /i = 0.260. Note that the values of the intermittency exponent /i for Il^^\e/ {{e))) 
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and for Ilg (^'■^■'/((e))) turn out to be the same. Other parameters are hsted in table [2] 
and table |3l 

The relative differences A„ = [U^^\e'^'^y{{e)))-U^^\e/{{6)))]/ui''\6/{{e))) for r/rj = 
6.57 and 65.7 are given, respectively, in figures H] (b) and (c), in which closed circles (full 
lines) represent experimental (theoretical) mean values of A„. The error bar in these 
figures is the standard deviation of 100 hidden (not appeared in the figures) data points 
for An which locate between the adjacent data points for A„ appeared in the figures. 
These figures show that the mean relative difference A„ in the region of central part 
of PDFs is about 10 times larger than the mean relative difference in the region of tail 
partH The small negative but nearly constant mean values of A„ in the tail part region 
tell us that the tail of Il^^\e / {{e))) and that of Il^^\e''^y {{e))) are parallel with each 
other, which gives the reason why we obtained the same value for both PDFs. We 
observe that the error bars in the tail part region become larger toward the rightmost 
end of PDF, which may be attributed to the smaller number of realizations in each bin 
there. Actually, the length of an error bar associated with a bin is quite close to the value 
^/T/WTT/N(^ which estimates the standard deviation of U^^\e^'^y {{e))) /u'i'^e / {{e))) 
with the help of the number of the realizations (A^*^*^)) in the bin under consideration 
for rig"'' (£:/((£:))) {Il^^\e^^'' / {{e)))) . The numbers of realizations in figure H] (c) are, for 
example, A^ = 21, A^(°) = 23 for the rightmost error bar, A^ = 206, A^(°) = 192 for 
the fifth error bar from the rightmost error bar, A^ = 31508, N^°^ = 31327 for an error 
bar at e/{{e)) ^ 5 and A^ = 2861811, A^(°) = 2886837 for an error bar at around the 
peak point of Il^^\e / {{e))) where £/{{£)) ~ 0.25. The e-dependence of the mean values 
of An in the central region indicates that the central part of Il^^\e^^'' / {{e))) around 
its peak point moves to rightwards relative to the central part of Il^\e / {{e))) . On 
the other hand, from the ^-dependence of the mean values of A„ in the tail region, it 
may be appropriate to interpret that the tail part of Ill^\e^'^y {{e))) moves to leftwards 
relative to the tail part of Il^^\e / {{e))) . If the number of realizations in each bin are 
increased, i.e., statistics are raised, we expect that the standard deviations of A„ should 
reduce their values and that the fluctuation of the mean values of A„ in the tail region 
should disappear. The difference of the squared time derivatives of ( |T3|) and ( |T2l) gives 
{A^^h/Atf - (Z\(3)t,/Zit)2 = {dv/dt){d'^v/dt'^)At + 0{Aty. From the direction of 
the relative horizontal shift of the PDFs, we know that the net contributions of the 
velocity component v for the region around the peak point and of the tail region satisfy, 
respectively, 

(dv/dt) {d\/dt^) > and (dv/dt) {d'^v/dt^) < 0. (14) 

Taking into account the smallness of the gradient of tail part PDFs, we see that the 
absolute value of the latter in (fT^ is quite large compared with the former value. 

The dependence of a* and 6 on r/rj {= in/v) ^-^^ given, respectively, in figure [3] 
(b) and (c) by open circles for 5 = 2, by open squares for 6 = 3 and by open triangles 

* Note that the connection pomts of H^"' (e/((e))) (H^"^ (e(°)/((e}))) for r/r] = 6.57 and 65.7 locate, 
respectively, at e*/{{e)) = 4.06 {e^°>/{{e)) = 4.17) and e*/((e)) = 5.20 ie^°>/{{e)) = 5.20). 
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for (5 = 5, which are extracted from the series of PDFs cleated with (I13p . The dashed 
hne in each figure, (b) and (c), represents an empirical formula obtained from all the 
data points for 6 = 2, 3 and 5 by the method of least squares. These figures prove 
again, even for the case of full contamination, the correctness of the assumption that 
the fundamental quantities of turbulence are independent of 5. We also found that ws is 
independent of 6 and has a common line log^gWs = 0.318 log^Q (r/r/) + log^g 0-141. The 
value of q' is found to be about q' = 1.05 (see table [2] (b)). 

There is only a slightly visible difference of the lines for a*, ws and of the values 
q' between those obtained from the two kinds of PDFs, one with less contamination 
( IT2l) and the other with full contamination (fT3l) (see figure [3] (b); see also table [2] and 
table [3]). The significant difference appears in the r/rj dependence of 6 which are shown 
in figure [3] (c). The difference in 6 explains the shift of the peak points between the two 
PDFs (see figures H (d)-(f)). 

5. Summary and Discussion 

The new scaling relation (|2]) is essential for the parameters ao, X and q, associated with 
the tail part PDF, to be determined self-consistently as functions of the intermittency 
exponent fi, and to be independent of the magnification rate S. On the other hand, 
we introduced the trial function ( ITTj) for the central part PDF with three adjustable 
parameters q', and 9, and found that these parameters are also independent of 6, 
and satisfy scaling behaviors in wider area not restricted to the inertial range. 

The independence of n from 6 ensures the uniqueness of the PDF of a for any value 
of 6. The comparison between the empirical formulae for n given in figure |3] (a) and the 
theoretical formula ([6]) provides us with the estimation £o = 20.6 cm which is about the 
same as the correlation length 17.9 cm or the separation 20 cm of the axes of adjacent 
rods forming the grid. Here, we are assuming that the empirical formulae are effective 
even for r/t] <2 (see the discussion in section |3] about the effective region oir/rf). Note 
that C.Q gives an estimation of the diameter of the largest eddy within the energy cascade 
model. 

As for the parameters appeared in the trial function for the central part PDFs, 
exp[— (7(,^„)] in f lTTj) . the discoveries that q' ~ 1.05 and that and In 1^3 reveal scaling 
properties are quite attractive for the research looking for the nature of the fluctuations 
surrounding the coherent turbulent motion of fluid. The fact that the value q' is 
quite close to 1 indicates that the HCT type function in ( iTTl) . i.e., the part giving 
exp[— (7(,^„)](,^*/,^„)^~^, is close to an exponential function. There is no theoretical 
prediction yet, which is based on an ensemble theoretical aspect or on a dynamical 
aspect starting with the N-S equation, to produce the formula for the central part 
PDF that represents the contributions both of the coherent turbulent motion providing 
intermittency and of incoherent fluctuations (background flow) around the coherent 
motion. If one could succeed to formulate a dynamical theory which produces properly 
the formula for the central part of PDFs starting with the N-S equation, it may provide 
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us with the physical meaning of the parameters g', 9 and W3, and with an appropriate 
pathway to the dynamical approach, e.g., the renormalization group approach, to fully 
developed turbulence. A study to this direction is in progress. 

Introducing two difference formulae f|T2|) and (1131) for the estimate of dv/dt, i.e., 
A^'^^v/ At with less contamination and A'^^'^v/ At with full contamination, we performed 
a trial for the extraction of information from PDFs by comparing two kinds of PDFs 
for energy dissipation rates, Iif\e/(ie))) and Iif\e^'^^ /is))) with e oc [A^'^^ V / At)'^ and 
e^^^ oc {A'^'^'> V / AtY . We observed that the intermittency exponents for the two kinds 
of PDFs turn out to take the same value /i = 0.260 (see table [2] and table [3] for other 
parameters). Through the accurate analyses of PDFs, it was also revealed that the 
parameters for Ii^^\e / lief)) and for Ii^.^\e'^^^ / lle'f)) are independent of 5 thanks to the 
new scaling relation ([2]), and that they show quite similar scaling behaviors extending to 
the regions with smaller and larger r/rj values outside the inertial range (see figure E]). 
The connection points a* of the tail and central parts of the PDFs take almost the same 
value for each r /t] (see tableland figure [3] (b)). It is found that, among the parameters 
controlling the central part, only 9 has a relatively larger deviation between the two 
different PDFs (see table |2] and figured (c)), which is related to the shift of the peak 
point occurred between the two kinds of PDFs. Other parameters q' and do not have 
significant difference among the two PDFs (see table |2]). 

Observing the relative difference A„ between I{^^\e^'^y {{ej}) and I{^^\e / {{ej)) in 
figures Kb) and (c) with the values </((e„)) = 4.53 for r/r] = 6.57 and e*J {{Sn)) = 5.20 
for r/r) = 65.7, we notice that the connection point e* of the center part PDF and the 
tail part PDF provides us with the boundary dividing two regions according to their 
nature of stability specified by the inequalities in ( !T4l) . It seems to tell us that the 
net behavior of incoherent motion of fluid contributing mainly around the peak point 
(central part) of PDF is an unstable time-evolution, whereas that of coherent turbulent 
motion contributing mainly to the tail part of PDF is a stable time-evolution. The 
former may be attributed to a manifestation of fluctuations, whereas the latter to the 
characteristics of intermittency. Note that we assumed that the central part n3"''j.(^n) is 
constituted of two contributions, one from the coherent contribution n3"g(£:„) and the 
other from the incoherent contribution AIl^^\en), and that almost all the contribution 
to the tail part fl^^^i{C,n) comes from the coherent intermittent motion of turbulence. 
Further investigation about these outcomes and their interpretation is necessary, which 
we leave as one of the attractive future problems]^ 

Let us close this paper by noting the studies in progress which are deeply related to 
the present work. It has been revealed (Motoike and Arimitsu 2012) that the new scaling 
relation ([2]) is intimately related to the 6°° periodic orbits {6 > 3) located at the threshold 
to chaos via the 5 (> 3) times ramification (bifurcation) in 5-period windows in the 

tt We observed that there is no visible difference between n3"''(e/((e))) and the PDF extracted with 
the formula dv/dt ~ [v{t + At) — v{t — At)]/2At which is correct without contamination up to the term 
oiO{At). 
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chaotic region, for example, of the logistic map. The self-similar nesting structure of 5^- 
period windows {k = 1, 2, 3, • ■ ■) can be an origin of intermittent coherent motion in fully 
developed turbulence. We expect that further investigation to this direction to extract 
a message provided by the new scaling relation may lead us to a novel interpretation of 
fully developed turbulence. We are also performing a precise comparison between the 
results extracted in this paper for the grid turbulence in a wind tunnel and those for 
4096^ DNS turbulence by raising the resolution of PDFs, i.e., by creating more data 
points for PDFs. The results of these studies will be published elsewhere in the near 
future. 
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